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A review is given on theoretical studies concerning the pairing symmetry in organic supercon- 
ductors. In particular, we focus on (TMTSF) 2 X and k-(BEDT-TTF) 2 X, in which the pairing 
symmetry has been extensively studied both experimentally and theoretically. Possibilities of 
various pairing symmetry candidates and their possible microscopic origin are discussed. Also 
some tests for determining the actual pairing symmtery are surveyed. 

KEYWORDS: organic superconductivity, pairing symmetry, (TMTSF) 2 X, k>(BEDT-TTF) 2 X, spin and 
charge fluctuations, Fermi surface 



1. Introduction 

Possible occurrence of unconventional superconductiv- 
ity in organic conductors 1-3 has been of great interest 
recently. Microscopically understanding the mechanism 
of pairing in those materials is an intriguing theoreti- 
cal challenge. Among the various candidates of uncon- 
ventional superconductors, in this paper we will focus 
on two groups of superconductors in which the pairing 
symmetry has been extensively studied both theoreti- 
cally and experimentally, namely, (I) (TMTSF) 2 X, 4 or 
the Bechgaard salts, where TMTSF is an abbreviation 
for tetramethyltetraselenafulvalene and X stands for an 
anion such as PF 6 , AsF 6 , C10 4 , etc., and (II) k-(BEDT- 
TTF) 2 X, 5,6 where BEDT-TTF is an abbreviation for 
bisethylenedithio-tetrathiafulvalene and X=Cu(NCS)2, 
Cu[N(CN) 2 ]Br, Cu 2 (CN) 3 , I 3 , etc. The key factors to 
be focused throughout the paper are the band struc- 
ture and the shape of the Fermi surface, the band fill- 
ing, and the wave number dependent pairing interactions 
mediated by spin and/or charge fluctuations and/or by 
phonons. Superconductivity near charge ordered state as 
in 6»-(BEDT-TTF) 2 X and a-(BEDT-TTF) 2 X 7 has also 
been investigated extensively, but will not be discussed 
here. 8 Superconducting states induced under high mag- 
netic fields, such as the Fulde-Ferrel-Larkin- Ovchin- 
nikov (FFLO) state, 9, 10 are also beyond the scope of the 
present paper. 11 

2. (TMTSF) 2 X 

2.1 Lattice Structure and the Phase Diagram 

The lattice structure of (TMTSF) 2 X is shown in Fig.l. 
The molecules are stacked along the a-axis (denoted as a 
hereafter), which is the most conducting axis because the 
overlap of the molecular orbitals, oriented in the stacking 
direction, is large. The molecules are weakly dimerized 
along the stacks. The charge transfer with the anions 
existing in between the conducting stacks results in one 
hole per two molecules. There is a weak overlap of the 
orbitals in the b direction, resulting in a weak two di- 
mensionality. 

A schematic phase diagram of (TMTSF) 2 X is shown 
in Fig.l. At ambient pressure, (TMTSF) 2 PF 6 undergoes 
a 2/c_f spin density wave (SDW) transition at 12K. Upon 
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Fig. 1. Left panebLattice structure of (TMTSF) 2 X in the o- 
b plane. Upper right panel: Typical shape of the Fermi sur- 
face. (Reprinted with permission from ref.24. Copyright 1983 by 
EDP Sciences.) Lower right panel: Schematic phase diagram of 
(TMTSF) 2 X. 

increasing hydrostatic pressure, the SDW transition tem- 
perature decreases, and superconductivity with a tran- 
sition temperature (T c ) of 0.9 K appears at 12 kbar. 12 
A similar phase diagram is obtained for X=AsF 6 . 13 It 
should be mentioned here that X-ray diffuse scattering 
experiments have revealed a coexistence of 2kp charge 
density wave (CDW) in the SDW phase for X=PF 6 , 14, 15 
while the amplitude of the 2kp CDW is very small for 
X=AsF 6 . 15 It should also be noted that the easy axis of 
the SDW is in the b' direction, 16, 17 which is the direction 
normal to the a-c plane and somewhat tilted from b due 
to the triclinic symmetry of the lattice. 

(TMTSF) 2 C104 becomes superconducting at ambi- 
ent pressure when the system is cooled down slowly 
enough for the anions to order at 24 K. 18 On the other 
hand, when the cooling rate is fast, the anions are 
frozen in random directions, and in this case, SDW takes 
place instead of superconductivity. 19 It has also been re- 
vealed that superconductivity is destroyed upon alloying 
(TMTSF) 2 C10 4 with a small amount of Re0 4 , and with 
further alloying, an SDW phase appears. 20-23 



2 J. Phys. Soc. Jpn. 

2.2 Electronic Structure 
Reflecting the lattice structure and also the anisotropy 

of the orbitals, the band structure of (TMTSF) 2 X is 
strongly one dimensional, i.e., the ratios of the hopping 
integrals in a, b, and c directions are tb/t a ~ 0.2 and 
t c /t b ~ 0.05, where t a = 200 - 300 mcV. 24,25 Since t c 
is extremely small, it is highly likely that the essential 
mechanism of the superconductivity lies within the two 
dimensional lattice {a-b plane). A typical Fermi surface 
is shown in Fig.l, which is open in the kb direction due 
to the quasi-one-dimcnsionality. We stress here that the 
anisotropy of the hopping integrals within the a-b plane 
largely owes to the fact that molecular orbitals are di- 
rected toward the a direction, while the distance between 
the molecules in the b direction is only about two times 
larger than that in the a direction. The hopping integral 
t a alternates along the a direction by about 10 ~ 20% 
due to the dimerization of the molecules. If we neglect 
this dimerization, the system is described by a 3/4- filled 
single band model, whose band dispersion is given as 

e(k) = 2t a cos(k a ) + 2tbCOs(kb). (1) 

Here, only the hoppings between the nearest neighbor- 
ing molecules in the a and b directions are considered. 
Lattice constants (neglecting the dimerization) are taken 
as the units of the length. Many of the theoretical ap- 
proaches have been based on this 3/4- filled band model, 
but in some studies, the strong dimerization limit has 
been assumed, where each dimcr of molecules is con- 
sidered as a site, so that the band now becomes half 
filled. 26 ' 27 There is also a study based on a two band 
model that maintains the realistic dimerization struc- 
ture. 28 

2. 3 Experiments Concerning the Pairing Symmetry and 
Their Theoretical Interpretations 

Early experiments for (TMTSF) 2 X, such as the spe- 
cific heat 29,30 and the upper critical field measurc- 
ments 31 ~ 33 had been interpreted within the conventional 
s-wave pairing. However, Abrikosov 34 pointed out the 
possibility of spin-triplet pairing based on the fact that 
T c is very sensitive to the existence of non-magnetic de- 
fects. 20, 21,35 ' 36 More recently, Joo et al. 22 * 23 have shown 
that the sensitivity of T c to non-magnetic impurities 
(Re0 4 ) in (TMTSF) 2 C10 4 is precisely what is expected 
from the T c reduction formula 37,38 for unconventional 
pairing. Generally, in a superconducting state that sat- 
isfies the condition 

JV(fe,iw n ) = 0, (2) 
k 

the presence of non-magnetic impurities is pair breaking, 
and thus strongly suppresses T c . 39 Here, F is the anoma- 
lous Green's function, and the condition (2) roughly cor- 
responds to a vanishing summation of the superconduct- 
ing gap function A(fe) over the Fermi surface. The T c 
reduction in this case is given in the form, 38 
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where ip is the digamma function, T c o is the transition 
temperature without impurities, and a is the pair break- 
ing parameter, which is determined by the scattering rate 
due to non-magnetic impurities. This is in fact the same 
as the formula derived by Abrikosov and Gor'kov for the 
case of s-wave pairing with magnetic impurities. 37 Since 
a triplet superconductivity has an odd parity gap, T c 
should be sensitive to the introduction of non-magnetic 
impurities. Note, however, that the condition (2) can be 
satisfied for a superconducting state with an even parity 
gap that changes sign on the Fermi surface, so that the 
sensitivity to the presence of impurities alone of course 
does not necessarily imply triplet pairing. In this sense, 
the sensitivity of the T c to non-magnetic defects concerns 
the orbital part of the pair wave function. 

Another experiment that indicated the possibility of 
unconventional pairing concerning the orbital part is the 
NMR experiment for X=C104 performed by Takigawa 
et al. 40 Namely, the proton spin-lattice relaxation rate 
l/Ti at zero magnetic field exhibits no coherence peak, 
and follows a power law temperature dependence close 
to T 3 . Such a behavior is generally characteristic to su- 
perconductivity with a gap having line nodes. 39 In fact, 
Hasegawa and Fukuyama 41 studied various types of sin- 
glet and triplet anisotropic pairings within the mean 
field approximation for a model with on-site and near- 
est neighbor attractive interactions, and calculated the 
spin-lattice relaxation rate. There it was shown that a 
singlet pairing without gap nodes on the Fermi surface 
exhibits a large coherence peak followed by an exponen- 
tial decay of l/Ti, while for a triplet pairing with gap 
nodes at k a = and thereby no nodes on the Fermi sur- 
face, which will be called p- wave hereafter (Fig. 2(c)), the 
coherence peak becomes smaller but still exists. Singlet 
and triplet pairings with line nodes of the gap intersect- 
ing the Fermi surface cannot be distinguished from the 
temperature dependence of l/Ti; they both exhibit es- 
sentially no (or very small) coherence peak and a power 
law decay roughly proportional to T 3 , which is similar 
to the experimentally observed behavior. 40 A more re- 
cent l/Ti measurement has been performed on X=PFg 
by Lee et a/., 42,43 who have found a similar behavior of 
l/Ti when a small magnetic field H is applied parallel 
to b', but also an anomalous l/Ti ~ T at low tempera- 
tures for high magnetic fields. On the other hand, Bclin 
and Behnia showed for X=C104 that the thermal con- 
ductivity rapidly decreases with lowering the tempera- 
ture below T c , indicating the absence of low lying excita- 
tions, and thus a fully gapped superconducting state. 44 
A possible explanation for this discrepancy between the 
conclusions of the NMR and the thermal conductivity 
experiments will be discussed in section 2.7. 

Before discussing the experimental results concerning 
the spin part of the pair wave function, let us briefly sum- 
marize some general aspects of spin triplet pairing. 39 In 
the case of triplet pairing, both the diagonal and the 
non-diagonal elements of the superconducting order pa- 
rameter matrix 

A(k)-( A " (fe) A ^ {k) ) (4) 
A iT (fc) A u (fc) ) (4) 
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remain finite in general, where A CT(T / are given as Aff = 
-d x + id y , A n = A iT = d z , and A u = d x + id y , using 



the vector d — (d x ,d y ,d z ). The order parameter vector 
d lies in the direction perpendicular to the total spin of 
the triplet pairs. 

Spin-triplet superconductivity can be identified by the 
NMR Knight shift measurement, which probes the uni- 
form spin susceptibility. The Knight shift decreases below 
T c for singlet pairing, while it stays constant for triplet 
pairing when the magnetic field H is applied perpendicu- 
lar to d. Another possible way of detecting triplet pairing 
is to measure the upper critical field H c2 . Cooper pair- 
ing under magnetic field is limited by both orbital and 
paramagnetic effects. 45 ' 46 For triplet pairing, however, 
the paramagnetic limit (the Pauli limit, or the Clogston- 
Chandrasckhar limit) is overcome when the magnetic 
field is applied perpendicular to d. 

Now, as for the actual experimental results, Lee et al. 
found for X=PF@ that the Knight shift does not de- 
crease below T c for magnetic fields applied parallel to 
a 42 (1.43T) or b' (2.38T). 43 These results indicate that 
the pairing indeed occurs in the spin-triplet channel, and 
that either d | c, or d rotates in accord with the direction 
of the magnetic field to satisfy d _L H. 

As for the upper critical field H c2 , Gor'kov and Jerome 
pointed out in the early days that H c2 extrapolated to 
T = may largely exceed the Pauli limit, suggesting the 
possibility of spin-triplet pairing. 47 More recently, the 
upper critical field has been studied with higher accu- 
racy and with precise orientation of the magnetic fields. 
H^ 2 , the upper critical field for H \\ b' has been found to 
exceed the Pauli limit for X=PF@ 48 and also for CIO4. 49 
Even if a spin-triplet pairing occurs, the pairing can still 
be orbitally limited, but Lebed 50 and later Dupuis et 
al. 51 showed that a magnetic field induced dimensional 
crossover from three to two dimensions can strongly en- 
hance the orbital limit of the critical field. Thus, as far 
as H h c2 is concerned, the experimental results seem to 
be consistent with the above interpretations of triplet 
pairing with d || c, or a rotatable d. 52 However, the in- 
terpretation on the temperature dependence of H " 2 has 
been controversial. For X=PF6, there is an inversion be- 
tween H® 2 and H^ 2 , where H" 2 > H h c2 for small magnetic 
field, but H« 2 < H b c2 for H > 1.6T. Moreover, ff c a 2 (T) 
as a function of T changes from a convex to a concave 
curve above H = 1.6T. From these experiments, Lebed 
et al. proposed that d b = and d a =^0 (d = (d a ,d b ,d c )) 
assuming strong spin-orbit coupling, so that the pairing 
is Pauli-paramagnetically limited for H \\ a for H < 1.5 
T, while the change of the curvature of H^ 2 {T) for higher 
magnetic fields may be because d rotates to become per- 
pendicular to H, or may be due to an occurrence of the 
FFLO state. 53 On the other hand, Duncan et al. argued 
that spin-orbit coupling should be weak since the heav- 
iest element in (TMTSF) 2 X is Se, so that d should be 
able to rotate according to the direction of H even for low 
magnetic fields. 54 A clear understanding for the direction 
of d, provided that triplet pairing does indeed occur, 55 
requires further theoretical and experimental study. 



2.4 Spin-fluctuation-mediated d-wave pairing 

In this and the next two subsections, we discuss some 
mechanisms for anisotropic pairing in TMTSF salts. 
Since the superconducting phase lies close to the SDW 
phase, and a number of experiments suggest the possi- 
bility of anisotropic, unconventional pairing, it is natu- 
ral to expect that the spin fluctuations mediate (or at 
least play an important role in) the Cooper pairing in 
TMTSF salts, as was pointed out by Emery. 56 The spin- 
fluctuation-mediated pairing scenario has in fact been 
supported by several theoretical studies on the quasi-one- 
dimensional Hubbard model, in which the on-site repul- 
sive interaction U is considered along with the kinetic 
energy part considered in section 2.2. The Hamiltonian 
is given in standard notation as, 

<i,j>,<7 i 

where c\ a creates an electron with spin a at site i (i.e., 
the i-th molecule), — c\ a Ci a , and tij — t a and tij = % 
for intrachain and interchain nearest neighbor hoppings, 
respectively. There, superconductivity has been studied 
using random phase approximation (RPA), 57 fluctua- 
tion exchange approximation (FLEX), 26 third order per- 
turbation, 27 or quantum Monte Carlo method. 58 Here, 
based on RPA equations (which will be written down in 
a general form for later use) for the single band Hub- 
bard model at quarter filling (quarter filling of holes, 
to be precise), we summarize the mechanism in which 
2k f spin fluctuations lead to ci-wave like pairing. Within 
RPA, singlet and triplet pairing interactions are given in 
the form, 59 - 61 

V s (q) = U + V(q) + ^U 2 x s (q)-±(U + 2V(q)) 2 x c (q) 



V\q) = V(q) 



\u 2 Xs {q)- l -(U + 2V(q)) 2 X c{q), (5) 



where V(q) is the Fourier transform of the off-site inter- 
actions (electron interactions between nearest neighbors, 
etc.), which is for the Hubbard model. Here, %s and Xc 
are the spin and the charge susceptibilities, respectively, 
which are given as 

Xo(<?) 



Xs(q) 



Xc(q) 



l-Uxo(q)' 

xo(q) 



l + (U + 2V(q)) X0 (q)' 
Here %o is the bare susceptibility given by 

1 ^f(e(p + q))-f(e(p)) 



(6) 



xo(q) 



p 



e(p) - e(p + q) 



where /(e) is the Fermi distribution function. Within the 
weak coupling BCS theory, T c is obtained by solving the 
linearized gap equation, 

, tanh(/3e(fe')/2) 



A 5 '*A 5 '*(fe) 



■£>«■«(*-*')- 
k' 



2e(k') 



A s >'(fc'). 
(7) 
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The eigenfunction A 5 '* of this eigenvalue equation is the 
gap function. The transition temperature T c is deter- 
mined as the temperature where the eigenvalue A reaches 
unity. In the summation over k' in the right hand side of 
eq.(7), the main contribution comes from k' on the Fermi 
surface because of the factor tanh (fie(fc )/ 2 ) jf we multi- 

2e(fe') 

ply both sides of eq.(7) by A s ' t (k) and take summation 
over the Fermi surface, we see that the quantity 

J2 kk/ v^(k-k')A^(k)A^(k') 

Eki^mr ( ' 

has to be positive and large in order to have large A, i.e., 
in order to have superconductivity with the gap A s ' t (k). 

Due to the good nesting of the Fermi surface, the 
bare susceptibility Xo{o) peaks at the nesting vector 
q = Q 2 fc F : an d since U > and V(q) = 0, Xs{q) becomes 
large at q = Q 2 fc F j while x c (q) remains small at all q. 
Within this formulation, the SDW transition tempera- 
ture is determined as the temperature where Uxo(Q2k F ) 
reaches unity. Thus, in the vicinity of the SDW transi- 
tion, the pairing interactions roughly satisfy the relation 
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(9) 



because the contribution from the spin fluctuations 
strongly dominates in eq.(5). Now, since the pairing in- 
teractions have large absolute values at q = Q 2 k F i the 
condition to have a positive V e g in eq.(8) can be approx- 
imately reduced to 

^ s '*(Q 2fcF )A s ' t (fe)A s '*(fe + Q 2kF ) < 0, k, k + Q 2kF e F.S. 



Fig. 2. Upper panel: candidates for the gap function of 
(TMTSF)2X are schematically shown along with the Fermi sur- 
face (solid curves), (a)d-wave, (b)/-wave, (c)p-wave. The dashed 
lines represent the nodes of the gap, whose kf, dependence is 
omitted for simplicity. "+,— " represent the sign of the gap func- 
tions. Lower panehFLEX calculation results for the two band 
model with finite dimerization. 28 (d)c(-wave gap, (e) /-wave gap, 
(f)|G(fc, i7rfcsT)| 2 , whose ridges represent the Fermi surface. 



(10) 

From this condition and eq.(9), we can see that the gap 
function has to change sign between k and k + Q 2 k F 
for singlet pairing, while the sign has to be the same 
across the nesting vector for triplet pairing. Since the 
spin part of the pair wave function is antisymmetric 
(symmetric) with respect to the exchange of electrons 
for spin singlet (triplet) pairing, the orbital part of the 
wave function, namely the gap function, has to satisfy 
the condition A s (k) = A s (— k) (even parity gap) and 
A t (k) — — A'(— ft) (odd parity), for singlet and triplet 
pairings, respectively. The gap functions satisfying these 
conditions are schematically shown in Fig. 2(a) (b). We 
will call the singlet pairing "ci-wave" in the sense that 
the gap changes sign as H 1 — along the Fermi sur- 
face, while the triplet pairing will be called "/-wave" in 

the sense that the gap changes sign as H 1 1 — . 62 

Since the pairing interaction is three times larger for the 
singlet pairing, c?-wave pairing takes place in this case. 
Note that a simpler form of an odd parity gap is the p- 
wave shown in Fig. 2(c), which changes sign as H — along 
the Fermi surface. However, this gap does not satisfy the 
condition (10) because the triplet pairing interaction is 
negative for the Hubbard model at least within RPA. 

Although we have adopted RPA equations in the 
above, similar conclusions have been drawn from other 
approaches as mentioned above. For example, an ap- 
proach along the line of RPA, but more suitable for deal- 
ing with strong spin fluctuations, is the FLEX method. 63 
In the FLEX, (i) Dyson's equation is solved to ob- 
tain the renormalized Green's function G(k), where k = 



(k, ie n ) denotes the wave vectors and the Matsubara 
frequencies, (ii) the effective electron-electron interac- 
tion V^\q) is calculated by collecting RPA-type dia- 
grams consisting of the renormalized Green's function, 
namely, by summing up powers of the irreducible sus- 
ceptibility Xirr(g) = — jj J2k + l)G(k) (A^number 
of fc-point meshes), (iii) the self energy is obtained as 
s O) = jf J2 q G(k - q)VW(q), which is substituted into 
Dyson's equation in (i) , and the self-consistent loops are 
repeated until convergence is attained. 

To obtain T c , the linearized Eliashbcrg equation for 
the singlet or the triplet gap function A s ' t (k), 

AA s '*(fc) = ~^2 V s ' l {k - fc')G(fc')G(-fc')AM(fc'), 

k' 

(11) 

is solved, where the singlet or the triplet pairing interac- 
tions V s,t are given again in the RPA form but using the 
irreducible susceptibility obtained from the renormalized 
Green's functions instead of the bare susceptibility. T c is 
the temperature where the eigenvalue A reaches unity. 

Kino and Kontani applied FLEX to the half-filled Hub- 
bard model, i.e., the model in the strong dimerization 
limit, 26 where they obtained a finite T c for the d-wave 
pairing. Kuroki et al. applied FLEX to a two-band model 
with finite dimerization and with next nearest neighbor 
interchain hoppings, and also found that the d-wave pair- 
ing (Fig. 2(d)) is the most dominant pairing, while triplet 
/-wave pairing (Fig. 2(e)) is subdominant. 28 

Another approach for the Hubbard model is the per- 
turbational theory, where all the Feynman diagrams up 
to a certain order are taken into account in the calcula- 
tion of the pairing interactions. Applying the third order 
perturbation theory to the half-filled model in the dimcr 
limit, Nomura and Yamada obtained finite values of T c 
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for the d-wave pairing. It has been found there also that 
the /-wave pairing is subdominant. 27 

As for numerical calculations for finite size systems, 
Kuroki and Aoki 58 adopted the ground state quantum 
Monte Carlo (QMC) technique. 64 ~ 66 This method en- 
ables us to accurately calculate correlation functions 
within statistical errors for finite size clusters. Since the 
superconducting order parameter is always zero for fi- 
nite size systems, we instead calculate its fluctuation, 
namely, the pairing correlation function, given in the 
form (ci + sCiCjC^ +s ) , where i,j denotes the sites, and i 
and i + 5 are the sites at which the Cooper pair is formed. 
When the tendency towards superconductivity is strong, 
the pairing correlation decays slowly at large distances 
between sites i and j. Applying this method to the sin- 
gle band Hubbard model at quarter filling, it has been 
found that the <i-wave pairing correlation function is en- 
hanced at large distances by the presence of the on-site 
repulsion U. 5& More recently, Kuroki et al. studied the 
pairing symmetry competition on the Hubbard model at 
quarter filling using the ground state QMC, where they 
found that d-wave and /-wave strongly dominate over 
p- wave. 

Apart from the studies directly dealing with the Hub- 
bard model, low energy theories using the interacting 
electron gas model like those for the purely one dimen- 
sional systems as will be mentioned in section 2.6 can be 
effective, but since the nodes of the d-wave gap run par- 
allel to the kt axis, it is necessary to take into account 
the quasi-one-dimensionality (the warping of the Fermi 
surface) to study d-wave pairing in a realistic situation. 
Duprat and Bourbonnais indeed showed the occurrence 
of d-wave pairing near the SDW phase within a renor- 
malization group study that takes into account the quasi 
one dimensionality. 68 

2.5 Spin triplet f-wave pairing 

Nevertheless, the spin-fluctuation-mediated d-wave 
pairing scenario contradicts with the experimental facts 
pointing towards spin-triplet pairing, especially for 
X=PFg. 42 ' 48 (Note that most of the d-wave theories ap- 
peared before the Knight shift measurements.) Kuroki et 
al. 28 provided a possible solution for this puzzle by re- 
calling that 2&f CDW actually coexists with 2k F SDW 
in the insulating phase for X=PF 6 . 14 ' 15 If 2k F CDW co- 
exists with SDW in the insulating phase, it is natural to 
assume that 2k f spin and 2k p charge fluctuations co- 
exist in the metallic phase lying nearby. Assuming the 
presence of charge fluctuations along with spin fluctua- 
tions with possible magnetic anisotropy (i.e., presence of 
easy and hard axes), the pairing interactions are given in 
generic forms, 

V s (q) = \v^{q)+V+-{q)-\v ch {q) 
V^{q) = -\v£(q)-\v*.{q) 

V'Hq) = \v: v z {q)-V+-{q)-\v ch { q ) (12) 

where V s z p z and are the contributions from lon- 

gitudinal and transverse spin fluctuations, respectively, 




with d//z 



Fig. 3. Phenomenological phase diagram of the pairing symmetry 
for (a)z || 6 28 and (b) z || c. 

while V c h is the contribution from the charge fluctua- 
tions. There are two triplet pairing interactions: V tl - for 
d _L z and V*" for d || z. The contribution from the spin 
fluctuations is expected to be large in the easy axis direc- 
tion of the SDW ordering. Then, taking the easy axis as 
the z-axis, we may assume V z p z (Q 2kF ) > V+-{Q 2kp ) be- 
cause the longitudinal spin susceptibility should exhibit 
stronger divergence at q = Q 2kp than the transverse 
ones near the SDW transition. Thus, —V {Q 2k ) > 

-V tll (Q2k F ) hol ds from eq.(12), where -V tJ -(Q 2kF ) is 
always positive. Furthermore from eq.(12), we can see 
that 

-V f± (Q 2kF ) > V s (Q 2kF ) (13) 

holds when the condition, 

V ch (Q 2kF ) > V+~(Q 2kF ) (14) 

is satisfied. This kind of relation between the sin- 
glet and the triplet pairing interactions when spin and 
charge fluctuations coexist has been pointed out by Taki- 
moto 69 for another candidate for a triplet superconduc- 
tor, Sr 2 Ru0 4 . 7 ° 

Now, an important point for a quasi-one-dimensional 
system is that the number of gap nodes that intersect 
the Fermi surface is the same between d- and /-waves 
due to the disconnectivity of the Fermi surface, so that 
which one of these two dominates is determined solely by 
the magnitude of the pairing interactions. Thus, triplet 
/-wave pairing with d perpendicular to the easy axis 
direction dominates over singlet d-wave when the contri- 
butions to the pairing interaction from the charge fluc- 
tuations is larger than that from the spin fluctuations in 
the hard axis direction. 

The above argument can be summed up as a phe- 
nomenological phase diagram shown in Fig. 3(a). In this 
phase diagram, there exists a region where p-wave pair- 
ing dominates because V*" {Q 2kF ) > holds when V z p z > 
2V^ - + V c h, namely, when the magnetic anisotropy is 
strong and the charge fluctuations are weak, so that a 
triplet gap that has different signs at both ends of the 
nesting vector can be favored. In this case, d is parallel to 
the easy axis. This p-wave mechanism has in fact been 
proposed for Sr2Ru04. 71, 72 s-wave pairing having the 
same gap sign over the entire Fermi surface is expected 
to dominate when the charge fluctuations are sufficiently 
strong because the singlet pairing interaction turns neg- 
ative (which is unrealistic for (TMTSF) 2 X). 
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(a) + d-wave + 

2k F SDw ^TOdJ^ O (D O © 

— p-wave + 
(with S z =0) 

(b) 

— f-wave + 
2k F CDW 

Fig. 4. (a)2fc F SDW configuration and (b) 2k F SDW+2fc F CDW 
configuration with likely pairings when the configuration "melts" 
to become metallic. 

An intuitive understanding for this phase diagram 
can be given as follows. In the 2k p SDW configura- 
tion, electrons (or, actually, holes in a 3/4- filled system) 
with antiparallel spins sit at next nearest neighbors as 
shown in Fig. 4(a), so if this configuration "melts" to be- 
come metallic, singlet pairing superconductivity with an 
even parity gap of A(fe) = + exp(i2k a ) + exp(—i2k a ) ~ 
cos(2/c a ) is likely to occur. "2" in the argument of "exp" 
implies that the pairs are formed at second nearest neigh- 
bors, and the "+" signs in front of the "exp" corre- 
sponds to singlet wave functions having the same sign 
in the right and the left directions, as shown in Fig. 4(a). 
Since the gap cos(2fc a ) has even parity and has nodes 
at k a — ±tt/4, this corresponds to the singlet d-wave. 
On the other hand, when 2k F SDW and 2k F CDW co- 
exist (namely when both SDW and CDW have a pe- 
riod of four lattice spacings), the electrons are aligned 
like in Fig. 4(b), so that when this configuration melts, 
triplet superconductivity with an odd parity gap of 
A(fe) = + exp(i4/c a ) — exp(— i4fc a ) ~ sin(4fc a ) is likely 
to take place. This corresponds to the /-wave gap. If we 
consider the magnetic anisotropy, a triplet pair formed 
at fourth nearest neighbors is expected to have a to- 
tal S z = ±1 because the SDW spins are oriented in 
the z direction, (z || b ), which explains d _L z for /- 
wave. On the other hand, if the pure 2kp SDW con- 
figuration (Fig. 4(a)) with z being the easy axis melts, 
a triplet pairing with S z = formed at next nearest 
neighbor sites may compete with the singlet pairing. 
This corresponds to the p-wave pairing with the gap 
A(fe) = +exp(i2fc a ) — exp(— i2k a ) ~ sin(2fc a ), whose 
nodes do not intersect the Fermi surface. Since the to- 
tal spin of a triplet pair in this case is expected to be 
perpendicular to z, d || z can be understood. 

In the above, the z-axis of the spins is taken in the 
b direction, namely, the easy axis direction, but if we 
assume that the spin fluctuations in the a-V plane are 
nearly isotropic and larger than those in the c (hard axis) 
direction, we can take the hard axis as the z-axis and 
thus V+-(Q 2kp ) > V s z p z {Q 2kp ), so that now d || c for /- 
wave pairing following a similar argument as before. This 
picture may be more suitable for (TMTSF^PFe since 
(i) the uniform susceptibilities for J? || a and H \\ b' are 
equal down to the very vicinity of the SDW transition, 
while that for H || c deviates from higher temperatures, 
and (ii) the direction of the SDW undergoes a spin-flop 
transition into the a direction under a magnetic field 




T/t a 

Fig. 5. Left: The model for (TMTSF) 2 X adopted in ref. 83. 
Right: The largest eigenvalue of the gap equation in the singlet 
and the triplet channels are plotted as functions of temperature 
for U = 1.7, V = 0.8, V = 0.45, V" = 0.2, V ± = 0.4, t b = 0.2, 
in units of t a . 83 

in the b' direction, 16, 17 which may be an indication that 
the 2kp spin fluctuations in the metallic state may not be 
so anisotropic within the a-b' plane. The phase diagram 
for z |j c assuming isotropic spin fluctuations in the a- 
b' plane is shown in Fig. 3(b). 73 Note that in this case, 
strong anisotropy in the spin fluctuations does not lead 
to p-wave pairing because the triplet pairing interactions 
always remain negative. 

After this phenomenological proposal and also a sim- 
ilar phenomenological argument of /-wave pairing by 
Fuseya et al., 7A studies based on microscopic models 
have followed. Tanaka and Kuroki considered a model 
which takes into account the off-site repulsive interac- 
tions X)<j j> Vij n i n j within the chains up to third near- 
est neighbors (Fig. 5, but with V± = 0), where the consid- 
eration of the second nearest neighbor repulsion V is the 
key. 75 This has been based on a consideration that since 
the coexistence of 2kp spin and 2kp charge fluctuations 
is necessary for /-wave pairing, and since the coexistence 
2kp SDW and CDW is experimentally observed, 14, 15 a 
model that can account for this coexistence should be 
the right Hamiltonian to be adopted. The mechanism 
of the coexistence of 2kp SDW and 2kp CDW itself had 
already been proposed by Kobayashi et al. and also stud- 
ied by Tomio and Suzumura, where the second nearest 
neighbor repulsion V' plays an essential role. 77-80 From 
Fig. 4, it can be seen how V induces the 2kp CDW in a 
quarter-filled system. When only the on-site U and the 
nearest neighbor V are present, the charges tend to take 
the 4fcp(= tt) CDW configuration, which has a period of 
two sites, while when V is present, the pairs of electrons 
sitting at second neighbors repel each other to result in 
the 2k F CDW(+SDW) configuration. 

For the Hamiltonian that considers U, V, V , and V" , 
the Fourier transform of the off-site repulsions, consid- 
ered in the RPA eq.(6), is given as 

V(q) = 2V cos(q x ) + 2V' cos(2 & ) + 2V" cos(3 fe ) (15) 

From eqs.(5),(6), and (15), it can be seen that 
Xs(Q 2 k F ) = Xc{Q 2 k F ) (where Q 2kp = (n/2,Tr)), and 
consequently — V t (Q 2 f. F ) — V s {Q2k F )' apart from the 
first order terms such as U + V(q), is satisfied when 
V = U/2. Within the phenomenological argument, this 
corresponds to the condition for /-wave to be degenerate 
with d-wave in the absence of magnetic anisotropy. 76 In 
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the actual RPA calculation for V' = U/2, /-wave slightly 
dominates over cZ-wave due to the effect of the first order 
terms neglected in the phcnomcnological argument. 75 

Fuseya and Suzumura 81 approached the same problem 
using the renormalization group method for quasi-one- 
dimensional systems along the line of Duprat and Bour- 
bonnais, 68 where a similar conclusion has been reached. 
Since the pairing competition is subtle, they further pro- 
posed a possible singlet d-wave to triplet /-wave transi- 
tion in the presence of magnetic field. 

According to the above studies, /-wave dominates 
over (i-wave when the second nearest neighbor repul- 
sion V' is equal to or larger than half the on-site re- 
pulsion U, which may be difficult to realize in actual 
materials. Nickel et al. have proposed a possible solu- 
tion for this difficulty, where they considered, in addi- 
tion to the intrachain repulsions, the interchain repul- 
sion and used the renormalization group technique for 
quasi-one-dimensional systems to reach a conclusion that 
/-wave dominates over (i-wave in a more realistic param- 
eter regime with a smaller second nearest neighbor re- 
pulsion. 4,82 Independently, Kuroki and Tanaka also con- 
sidered a model that considers the nearest neighbor in- 
terchain repulsion V± as shown in Fig. 5. 83 Within RPA, 
the term 2V± cos(q y ) is added in the right hand side of 
eq.(15), so that the condition for Xs{Q 2 k F ) = Xc{Q 2 k F ) 
now becomes V + V± = U/2. This is a much more re- 
alistic condition than V = U/2 because the interchain 
distance is similar to the intrachain second nearest neigh- 
bor distance, so that we can expect V± to be as large as 
V. The actual calculation shows that /-wave dominates 
(has a larger eigenvalue A) over d-wave for a parame- 
ter set, e.g., U = 1.7, V = 0.8, V = 0.45, V" = 0.2, 
V± = 0.4, tf, = 0.2 in units of i (Fig.5), where the rela- 
tive magnitude of the interactions can be considered as 
realistic. 

2.6 Other Mechanisms for Triplet Pairing: Phonons, 
Ring Exchange 
In this subsection, we discuss some other mechanisms 
for spin-triplet pairing proposed for (TMTSF) 2 X. From 
the early days, possibility of spin-triplet superconduc- 
tivity in (TMTSF)2X has been discussed in terms of 
the low energy effective theory called the g-ology ap- 
proach for the purely one dimensional interacting elec- 
tron gas, i.e., the Tomonaga-Luttinger model. 84 In the 
g-ology phase diagram, the spin-triplet superconducting 
phase and the SDW phase share boundary, 85 so that it is 
tempting to relate this superconducting state with that 
of (TMTSF) 2 X, as was discussed in some studies. 86 More 
recently, this phase boundary between the SDW and the 
triplet superconductivity has been discussed as having 
SO (4) symmetry. 87, 88 Nevertheless, since exact numeri- 
cal studies on the purely one-dimensional extended Hub- 
bard model, where the on-site U and the nearest neigh- 
bor V is considered, 89 show that superconductivity does 
not occur in a realistic parameter regime when the in- 
teractions are all repulsive, it is likely that some kind 
of attractive interaction, most probably originating from 
electron-phonon interaction, should be necessary in order 
to realize the triplet superconducting state in the g-ology 



phase diagram, as discussed in some studies. 90 ' 91 

Apart from the g-ology-type approach, there have been 
studies on the phonon mechanism of triplet pairing. 
Kohmoto and Sato proposed a p-wave pairing mechanism 
due to a combination of electron-phonon interaction, 21cf 
spin fluctuations, and the disconnected Fermi surface. 92 
Assuming that the electron-phonon interaction is weakly 
screened, a long-ranged attractive interaction arises in 
real space, which means that the pairing interaction be- 
comes large and negative around q ~ in momentum 
space. If we denote this interaction as —V e i- P h(q), and 
if the spin fluctuations also contribute to the pairing to 
some extent, the pairing interactions are given as 

V s (q) = -V cl - ph (q) + ^V s z ;(q) + V+-(q) 
VHQ) = -Ki-ph(<?) - \vg{q) 

V*Hq) = -V el - ph (q)+ 1 -V s l*(q)-V+-(q).(W) 

In ref.92, the competition between s- and p-wave pair- 
ings was discussed, while the possibility of d- and /-wave 
was not considered because the warping of the Fermi 
surface was neglected. Let us first neglect the magnetic 
anisotropy, i.e., V zz = . Around q ~ 0, neglecting 
the spin fluctuation contribution, the pairing interaction 
is negative and has the same magnitude between singlet 
and triplet pairings. Thus, s- and p-wave pairings, whose 
gap does not change sign on each portion of the discon- 
nected Fermi surface, are equally favored by this inter- 
action around q ~ 0. At q — Q 2 k F on the other hand, 
neglecting the electron-phonon interaction this time, the 
positive spin fluctuation contribution for the singlet pair- 
ing works destructively against s-wave because the gap 
does not change sign across Q 2/ t F , while the negative 
contribution for the triplet channel also works against 
p-wave, whose gap changes sign. Since this destructive 
spin fluctuation contribution is smaller for triplet pair- 
ing, p-wave dominates over s-wave. Note that here again, 
the close competition between p-wave and s-wave arises 
from the disconnectivity of the Fermi surface owing to 
the (quasi) one dimensionality, i.e, the additional node 
in the p-wave gap as compared to the s-wave does not 
intersect the Fermi surface. 

If we further take into account the magnetic anisotropy 
and assume V zz < by taking the hard axis (c-axis) 
in the z direction, the negative spin fluctuation contri- 
bution in the triplet pairing interaction is smaller (and 
thus favorable for p-wave pairing) for d _L z than for 
d || z. Therefore, if the direction of d for p-wave pairing 
is governed by the magnetic anisotropy of the SDW, d 
is likely to be perpendicular to the hard axis direction, 
namely, in the a-b plane for (TMTSF) 2 PF 6 . 

Suginishi and Shimahara also proposed a phonon- 
mediated mechanism for p-wave pairing. 93 By consider- 
ing moderately screened phonons and also including the 
corrections due to charge fluctuations, they obtained an 
attractive pairing interaction that has a large magnitude 
around q ~ and a small one around q ~ Q 2 fe F • By fur- 
ther considering the Coulomb pseudo potential, which 
suppresses only the s-wave pairing, it has been found 
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there that p-wave pairing dominates in a certain param- 
eter regime. 

Recently, Ohta et al. proposed a non-electron-phonon 
mechanism for spin-triplet pairing. 94 The mechanism is 
based on the fact that in a triangle lattice consisting of 
three sites with two electrons, a ferromagnetic interac- 
tion arises by considering a consecutive exchange of the 
positions of the electrons. 95 If ferromagnetic spin fluc- 
tuations arise due to this "ring exchange mechanism" 
on a certain lattice, triplet pairing superconductivity 
may take place. They considered the Hubbard model 
on a one-dimensional "railway-trestle" (or zigzag) lat- 
tice, where they used the density matrix renormaliza- 
tion group method to find that triplet pairing correla- 
tion functions decay more slowly than the singlet ones. 
Their numerical calculation has been restricted to purely 
one dimensional systems so far, but they further propose 
that this mechanism may be applicable to the quasi- 
one-dimensional material (TMTSF) 2 X since the signs of 
the intrachain and interchain hopping integrals (by con- 
sidering also the next nearest neighbor interchain hop- 
ping) 24, 25 satisfy the condition for the ferromagnetic in- 
teraction. 95 

In all the mechanisms discussed in this subsection, at 
least one of the 2fcp fluctuations, spin or charge, are not 
taken into account, although they should both be present 
at least for X=PF 6 . Then, whether both of these fluctu- 
ations play essential roles or not in the occurrence of 
superconductivity is the key toward clarifying whether 
/-wave discussed in section 2.5 or other triplet pairings 
dominate, provided that the pairing indeed occurs in the 
triplet channel. 

2. 7 Tests for the Pairing Symmetry Candidates 

In this section, we discuss some experimental tests (al- 
ready existing ones as well as proposals for future study) 
for the candidates for the pairing symmetry discussed 
above. For the pairing symmetries whose gap has line 
nodes on the Fermi surface such as d-wave and /-wave, 
the spin-lattice relation rate 1/Ti exhibits essentially no 
(or very small) coherence peak and a power law decay 
proportional to ~ T 3 , 96 which is consistent with the ex- 
periments for X=C104 40 and for X=PFg 43 at low mag- 
netic fields. On the other hand, whether these pairings 
can account for the peculiar behavior of 1 jT\ observed for 
X=PF6 at high magnetic field, i.e., 1/Ti ~ T at low tem- 
peratures as well as a small peak below T c for H | a, 42,43 
remains open as an interesting future study. 

At first glance, only p-wave and s-wave pairings seem 
to be consistent with the thermal conductivity measure- 
ment for X=C1C>4 suggesting a fully gapped state. 44 
However, Shimahara has argued that a fully gapped 
state is possible even for d-wave pairing particularly 
in (TMTSF) 2 C104 , because in this case, anion order- 
ing takes place above the superconducting T c , so that 
a "gap" opens up on the Fermi surface at positions 
(k a = ±7r/4) where the nodes of the superconducting 
gap would otherwise intersect (see Fig. 2(a)). 97 Exactly 
the same argument holds for /-wave pairing since the po- 
sitions of the gap nodes on the Fermi surface are the same 
between / and d. A fully gapped state usually results in 



a coherence peak followed by an exponential decay in 
1/Ti as mentioned in section 2. 3, 41 which seems to be 
in contradiction with refs.40 and 43, but since 1/Ti can 
be affected by the presence of impurities, 98 vortices, 99 
or correlation effects, 100 the clarification of the relation 
between 1/Ti and the thermal conductivity experiments 
is open for future study. 

From the microscopic view discussed in the preceding 
sections, /-wave and p-wave are the main candidates for 
spin triplet pairing. 42,43,48 As for the direction of d, if 
we assume that the spin fluctuations in the c direction 
are weak while those in the a-b planes have similar mag- 
nitude, d of /-wave pairing should lie in the c direction 
as discussed in section 2.5, which is consistent with the 
Knight shift results. 42, 48 If we assume on the other hand 
that the spin fluctuations are solely strong in the b' di- 
rection compared to those in the a-c plane, then /-wave's 
d is perpendicular to b' and lies in the a-c plane as also 
discussed in section 2.5, which is more closer to the d di- 
rection proposed by Lebed et al. 53 from the temperature 
dependence of and H^ 2 . 4S In the case of p-wave pair- 
ing, if the anisotropic spin fluctuations contribute to the 
pairing interaction in the form given in eqs.(16) (V e \- P h 
need not be due to phonons), d is likely to lie in the a- 
b plane. Thus, it may be possible to distinguish / and 
p from the direction of d, provided that the anisotropic 
spin fluctuations play a role in the Cooper pairing. Such 
a test, however, has to be done in the absence of, or un- 
der low, magnetic field since d may rotate regardless of 
the pairing symmetry if the magnetic field is sufficiently 
large to overcome the effect of the magnetic anisotropy. 

Although there exist few experiments up to date, pos- 
sibility of determining the pairing symmetry from tunnel- 
ing spectroscopy measurements has been proposed the- 
oretically by several groups. Sengupta et al. pointed out 
that the presence/absence of zero energy peak in the tun- 
neling conductance can be used to distinguish various 
types of pairings in (TMTSF) 2 X. 101 In fact, the zero en- 
ergy peak in the tunneling spectroscopies of anisotropic 
superconductors (those with sign change in the gap) orig- 
inates from the zero-energy Andreev bound state caused 
by the sign change of the pair potential felt by the quasi- 
particle in the reflection process at the surface, 1 2,10 
and has turned out to be a powerful method for prob- 
ing the pairing symmetry in anisotropic superconduc- 
tors such as the high T c cuprates. 104 In the case of tun- 
neling parallel to a in particular, the zero energy peak 
does not exist for d-wave. This is because the injected 
and the reflected quasiparticles feel the same gap due 
to A d (k a ,k b ) = A d (-k a ,k b ) (see Fig. 2(b)). By con- 
trast, the zero energy peak does exist for p-wave and 
/-wave, where Af^ p {k a , fc&) = —Af^ p (—k a ,kb) is satis- 
fied (Fig. 2(c)). Tanuma et al. further pointed out that 
p-wave and /-wave can be distinguished from the over- 
all shape of the surface density of states (overall struc- 
ture of the tunneling spectrum) because p-wave is a fully 
gapped state, which results in a U-shaped surface den- 
sity of states around the Fermi level, while /-wave results 
in a V-shaped one. 105 Therefore, the combination of the 
absence/presence of the zero energy peak and the overall 
shape of the spectrum enables us to distinguish p, d, and 
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/-wave pairings. 

Further theoretical studies based on various shapes of 
the Fermi surface have been performed. 106, 107 Tanuma 
et al. showed that when the Fermi surface is warped in 
a certain manner, the zero energy peak can appear even 
in the case of d-wave. In this case, d and /-wave can be 
distinguished by the way the zero energy peak splits in 
the presence of a magnetic field. 106 Such studies show 
that the existence of the zero energy peak is sensitive to 
the shape of the Fermi surface (compare Fig.2(a)and (d), 
or (b) and (e)). Since the hopping integrals, and thus the 
Fermi surface, of (TMTSF^X depend on the pressure, 
the temperature, and the anions, 24, 25, 108 it is necessary 
to strictly pin down the actual shape of the Fermi surface 
at the temperature and the pressure at which supercon- 
ductivity takes place in order to distinguish the pairing 
symmetry from the presence/absence of the zero energy 
peak. 

The tunneling tests above mainly concern the orbital 
part of the pairing. On the other hand, Bolech and Gi- 
amarchi proposed a tunneling experiment to distinguish 
directly the spin part of the pairing. 109 They showed that 
the I-V characteristics of a normal metal-triplet super- 
conductor junction are unaffected by an application of 
magnetic field perpendicular to d, while the Zeeman ef- 
fect affects the I-V characteristics when d \\ H similarly 
to the case of normal metal-singlet superconductor junc- 
tion. Therefore, the spin part of the pairing, whether it 
is singlet or triplet and also the direction of d if triplet, 
can be determined by measuring the I-V characteris- 
tics of the junction under a rotating magnetic field, pro- 
vided that d does not rotate according to the direction of 
the magnetic field. Vaccarella et al. also proposed a way 
of directly probing the spin part of the triplet pairing. 
Namely, they showed that the Josephson effect between 
two triplet superconductors is very sensitive to the di- 
rection of d across the junction, and proposed that this 
sensitivity can be used as a test for triplet superconduc- 
tivity. 110 

As a final remark in this subsection, it is important 
to recognize that the pairing symmetry might be dif- 
ferent for different anions. This possibility is suggested 
especially from the viewpoint discussed in section 2.5. 
For instance, the amplitude of 2kp charge fluctuations, 
which has to be large for /-wave to dominate over d- 
wave, is found to be small for X=AsFg, so that /-wave 
has less chance of dominating over d-wave than in the 
case of X=PF 6 . Thus, the pairing symmetry of a TMTSF 
superconductor with a certain anion should be deter- 
mined by a combination of multiple experiments on the 
salt with that very anion. Furthermore, we must keep 
in mind that the pairing symmetry might even change 
for the same salt under different environment, such as 
the pressure and the strength of the magnetic field, be- 
cause several pairing symmetries may be closely compet- 
ing. In particular, as mentioned in section 2.5, singlet 
to triplet transition may take place under high magnetic 
field since the singlet pairing is Pauli-paramagnetically 
limited. 81,111,112 
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Fig. 6. (a) The lattice structure of k-(BEDT-TTF) 2 X in the b-c 
plane. &i, & 2 , • " represent the hopping integrals in the four band 
model, (b) The lattice structure of the dimer model, (c) Phase di- 
agram of k-(BEDT-TTF) 2 X. 5 (d) Band structure and the Fermi 
surface of k-(BEDT-TTF) 2 Cu(NCS) 2 . (Reprinted with permis- 
sion from rcf.119. Copyright 1988 by the American Physical So- 
ciety.) 



3. k-(BEDT-TTF) 2 X 

3.1 lattice Structure and the Phase Diagram 

The lattice structure of k-(BEDT-TTF) 2 X is shown 
in Fig. 6 (a), which consists of dimers formed by a pair 
of face-to-face molecules, b- and c-axis are taken as in 
Fig. 6, while the BEDT-TTF layers and the anion layers 
alternate along the a-axis. Relatively large overlap of the 
orbitals between the dimers exists (see section 3.2) while 
the overlap between the BEDT-TTF layers is very small, 
resulting in a strong two dimensionality. 

In Fig. 6(c), the generic phase diagram of «-(BEDT- 
TTF) 2 X is shown, which has been extensively stud- 
ied by Kanoda et al. 5 The superconducting and the 
antiferromagnetic insulating phases share a first or- 
der phase boundary. Recently, this boundary has 
been revealed to persist above the superconduct- 
ing T c and the Neel temperature into the bound- 
ary of the paramagnetic insulating and the metallic 
phases, ending up at a certain critical point, 113, 114 
where an anomalous criticality has been found re- 
cently. 115 The horizontal axis in the phase diagram 
can be considered as hydrostatic or chemical pressure, 
where superconductivity with T c exceeding 10K occurs 
at ambient pressure for k-(BEDT-TTF) 2 Cu(NCS) 2 116 
and K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br, while k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]C1 is an antiferromagnetic insulator 
below 26K at ambient pressure 117 and becomes super- 
conducting with T c = 12. 8K under an applied pressure 
of 0.3kbar. 118 
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Fig. 7. (a) The d x i _ y 2-wave gap in the original Brillouin zone 
and in the unfolded one (right). Note that although the gap 
changes sign at the Brillouin zone edge (k c = ±7r), the nodes 
of the gap are not located there; the gap jumps from a positive 
to a negative value, as can be seen more clearly in the unfolded 
Brillouin zone, (b) The d xy -wave gap. Here, we show the case 
when the two portions of the Fermi surface splits due to the lack 
of center-of-inversion symmetry. In this case, the d xy gap nodes 
do not intersect the Fermi surface, although the gap does be- 
come small near the Brillouin zone edge. On the other hand, if 
the two portions stick, the d xy nodes intersect the Fermi surface 
at the Brillouin zone edge. The solid (dashed) curves represent 
the portions of the Fermi surface where the gap has a positive 
(negative) sign. Q represents the wave vector of the spin fluctu- 
ation mode that favors each pairing symmetry. The short arrows 
denote the positions of the gap nodes. 

3.2 Electronic structure 

The values of the intermolecular hopping integrals 
shown in Fig. 6 (a) have been estimated using the ex- 
tended Huckel method, 119, 120 which is summarized in 
Table I. The hopping integral in the 6-direction alter- 
nates as tbi, tb2, tbi, ■ ■ ■ , where \tbi\ > | it,2 1 because of 
the dimerization of the molecules. In Fig. 6(d), the band 
structure of k-(BEDT-TTF) 2 Cu(NCS) 2 119 is shown. 
Four bands exist near the Fermi level because there are 
four BEDT-TTF molecules per unit cell. Due to the 
dimerization, a gap opens up between the bonding and 
the antibonding bands. There is one hole per dimer, so 
that only the upper two bands cross the Fermi level, re- 
sulting in two portions of the Fermi surface (Fig. 7). The 
two portions of the Fermi surface are connected at the 
Brillouin zone edge for X=Cu[N(CN) 2 ]Br, Cu 2 (CN) 3 , 
±3, etc., in which the anions are arranged in a manner 
that the system possesses center-of-inversion symmetry, 
which results in t p = t p , t q = t' q . On the other hand, for 
X=Cu(NCS) 2 , the system lacks the symmetry so that 
tp 7^ tpi tq 7^ t'„, and in that case, the two portions of the 
Fermi surface are separated into an open Fermi surface 
and a closed one. 

In the limit of large tbi, namely, when the dimeriza- 
tion is strong, each dimer can be considered as a single 
site, so the system reduces to a two band model shown 
in Fig. 6(b) with n = 1, where the band filling is now 
defined as n=(the number of electrons/the number of 
sites). 121 > 122 In other words, the energy gap between the 
upper two and the lower two bands becomes large when 
the dimerization is strong, so that the lower two bands, 
which do not cross the Fermi level, can be neglected. In 
this strong dimerization limit, the effective hopping inte- 
grals U and t c are given as tb = — £fc 2 /2, t c = (—t p +t q )/2, 
and t' c = (— t' p + t' q )/2 in terms of the original hopping in- 
tegrals, 121 which gives \t b /t c \ ~ 0.8 for X=Cu(NCS) 2 and 



anion b± 62 p p' q g' 

Cu 2 (CN) 3 22.36 11.54 Koi - -2.90 ^~ 

Cu(NCS) 2 22.95 11.31 9.85 10.09 -3.30 -3.76 

Cu[N(CN) 2 ]Br 24.37 9.16 10.14 -3.40 

Table I. Hopping integrals estimated in ref.120. In units of 
10~ 2 eV. 

\t b /t c \ ~ 0.7 for X=Cu[N(CN) 2 ]Br. The system further 
reduces to a half- filled single band model when t c = t' c . 

3. 3 Experimental Results Concerning the Pairing Sym- 
metry 

Here, we summarize the experimental results concern- 
ing the pairing symmetry. 6 In the NMR experiments 
for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br, the 13 C Knight shift 
has been found to decrease below T c , 123, 124 which is con- 
sistent with singlet pairing. Also, the 13 C spin-lattice 
relaxation rate l/Ti for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br 
exhibits no coherence peak, and a power law decay pro- 
portional to T is seen below T c . 123 - 125 As in the case of 
(TMTSF) 2 X, this is consistent with the presence of line 
nodes in the superconducting gap. Also, in a thermal 
conductivity measurement for X=Cu(NCS) 2 , a T-linear 
term has been found at low temperatures, suggesting the 
existence of nodes in the gap. 126 

On the other hand, there has been much contro- 
versy concerning the measurements of other quantities. 
The magnetic penetration depth has been measured us- 
ing techniques such as muon spin relaxation, 127 ' 128 ac 
susceptibility, 129,130 surface impedance, 131 ' 132 and dc 
magnetization. 133, 134 The penetration depth should ex- 
hibit an exponentially decaying behavior for a fully 
gapped state, while a power-law dependence is expected 
at low temperatures for gaps with nodes. Some stud- 
ies have found for X=Cu(NCS) 2 that the tempera- 
ture dependence of the penetration depth is consis- 
tent with a conventional full gap state, 127 ' 131 ' 133 while 
others have found results consistent with a gap with 
nodes. 128 ' 129,132 ' 135 Similar controversy on the penetra- 
tion depth also exists for X=Cu[N(CN) 2 ]Br, where the 
presence of gap nodes 128 ' 130, 135 as well as the absence of 
them 131,133,134 has been suggested. 

The temperature dependence of the specific heat has 
also been another issue of controversy. Nakazawa and 
Kanoda found for X=Cu[N(CN) 2 ]Br a T 2 dependence 
of the electronic specific heat, 136 which was taken as an 
indication for the presence of nodes in the gap. How- 
ever, more recent results for X=Cu[N(CN) 2 ]Br 137 and 
for X=Cu(NCS) 2 138 have shown exponentially activated 
temperature dependence, indicating a fully gapped su- 
perconducting state. 

The above experiments do not give direct information 
on the position of, if any, the nodes in the gap function. 
Several groups have in fact made attempts to directly 
determine the node positions. A millimeter- wave trans- 
mission experiment suggested a gap function which has 
nodes in the direction shown in Fig. 7(a). 139 If we un- 
fold the Brillouin zone (right panel of Fig. 7(a)), which 
corresponds to adopting a single dimer as a unit cell, 
(this is possible when t c — t' c ), this gap function has 
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nodes in the diagonal direction like the gap in the high 
T c cuprates. 39, 104 In this sense, we will refer to this pair- 
ing symmetry as <i a ,2_ y 2-like pairing hereafter. (Note that 
this terminology is the opposite to the one adopted in 
rcf.140. Namely, we define the x and y axes by rotat- 
ing b and c axes by 45 degrees.) However, different in- 
terpretations on this experiment have been proposed af- 
terwards. 141 ' 142 On the other hand, Arai et al. showed 
for k-(BEDT-TTF) 2 Cu(NCS) 2 using in-plane STM mea- 
surement that the gap is the largest in the b and c direc- 
tions. 143 This is more consistent with the gap function 
shown in Fig. 7(b), which has nodes in between the kb and 
k c directions. Izawa et al. measured the thermal conduc- 
tivity of k-(BEDT-TTF) 2 Cu(NCS) 2 under a magnetic 
field rotating in the b-c plane, where an oscillation again 
consistent with the gap in Fig. 7 (b) was observed. 144 This 
conclusion is based on theoretical studies 145 showing that 
the density of states of a superconducting state having 
gap nodes in some direction becomes large (small) when 
the magnetic field is applied in the antinodal (nodal) 
direction due to the Doppler shift of the quasiparticle 
energy spectrum. 146 Since the gap in Fig. 7(b) has nodes 
in the vertical and horizontal directions in the unfolded 
Brillouin zone, we will call this gap d xy -\ike hereafter. 

3.4 Phonon Mechanisms 

If the symmetry is s-wave, the pairing is most likely 
due to electron-phonon interactions. Here we give a brief 
survey on the electron-phonon-interaction mechanism 
proposed for k-(BEDT-TTF) 2 X. Yamaji argued that in- 
tramolecular phonons (molecular vibrations) should play 
an important role in the occurrence of superconductivity 
in organic materials. The theory in which the coupling 
between the electrons and the molecular vibrations is 
taken into account was applied to /3-(BEDT-TTF) 2 I 3 147 
and to k-(BEDT-TTF) 2 Cu(NCS) 2 , 1 where T c for s-wave 
superconductivity was estimated. 

Girlando et al. 148 considered the coupling between 
the electrons and both the intramolecular and the in- 
termolecular phonons (lattice vibrations), and estimated 
the s-wave T c for k-(BEDT-TTF) 2 I 3 and /3-(BEDT- 
TTF) 2 I 3 using the Allen-Dynes formula. 149 They con- 
cluded that the contributions of both phonons are im- 
portant in understanding the experimental values of the 
T 

On the other hand, Varelogiannis considered the 
phonons with small q due to weak screening, 150 similar 
to those discussed in section 2.6 for TMTSF salts. Adopt- 
ing the single band dimer model, the BCS gap equation 
was solved by assuming the phonon-mediated attractive 
interaction in the form oc — l/(q 2 + q 2 ), where q c is a mo- 
mentum cutoff parameter, and also taking into account 
the Coulomb pseudopotential. There, it has been found 
that a close competition between an anisotropic s-wave 
and d x 2_ y 2-wa,ve pairings takes place. 

3.5 Studies on the dimer model 

If the superconducting gap in k-(BEDT-TTF) 2 X in- 
deed has nodes, the most probable scenario is that elec- 
tron correlation plays an important role in the pairing. In 



this and the next subsection, we discuss electronic mech- 
anisms of superconductivity having gap with nodes. 

Following the study by Kino and Fukuyama, who 
showed that the insulating state in k-(BEDT-TTF) 2 X 
can be understood as a Mott insulator of the half-filled 
dimer model in which the on-site repulsion U (which 
corresponds to the effective repulsion within the dimer) 
is considered, 5, 122 various theoretical studies concerning 
the superconductivity have been performed for the dimer 
model with on-site U on the two band 151 (for t c 7^ t' c ) 
or the single band (for t c = t' c ) lattice 58, 152-157 with 
tb/t c = 0.6 ~ 0.8. In these studies, the d x 2_ y 2-wa,ve has 
been found to be the most dominant pairing. Namely, 
in RPA or FLEX studies, 151-155 the spin susceptibility 
is found to have a peak at the nesting vector Q (near 
(it, it) in the unfolded Brillouin zone, see Fig. 7(a)) that 
bridges the open portions of the Fermi surface, although 
the nesting is not so good. Then, in order to have op- 
posite signs of the gap across Q, and also to satisfy the 
even parity condition for spin-singlet pairing, the d x 2_ y 2 
gap is favored. Within the FLEX studies, T c has been 
estimated to be O(10K), 152 ' 153 consistent with the ex- 
periments. By using the third order perturbation theory, 
Jujo et al. showed that d x 2_ y 2-w&ve pairing dominates 
but with T c lower than those obtained in FLEX. 156 They 
concluded that the vertex corrections that are not taken 
into account in FLEX have an effect of suppressing the 
T c especially for systems on frustrated lattices, at least 
up to third order. As for numerical approaches for finite 
size systems, Kuroki and Aoki applied the ground state 
QMC technique and showed that the d x 2 _ y 2 -wave pair- 
ing correlation function is enhanced at large distances ac- 
companied by a development of the spin correlation near 
q = (7r,7r). 58 Quite recently, Liu et al. used the varia- 
tional Monte Carlo technique, where d x 2_ y 2-w&ve super- 
conducting order parameter is found to be enhanced in 
a certain range of U/t. 157 

There have also been some strong coupling approaches 
along the line of Anderson's resonating valence bond 
(RVB) theory for the high T c cuprates, 158 in which large 
U/t is assumed. In the large U/t limit, the Hubbard 
model on a square lattice is transformed into the t-J 
model which consists of the nearest neighbor antifcrro- 
magnetic superexchange term (the magnitude of the su- 
perexchange being J = At 2 /U, where t is the nearest 
neighbor hopping) and the hopping term in the space 
that prohibits double occupancy of electrons at a sin- 
gle site. At exactly half filling, the hopping term van- 
ishes, so that the model reduces to the antiferromagnctic 
Heisenberg model, which describes the experimental sit- 
uation for the cuprates to be an antiferromagnetic Mott 
insulator when carriers are not doped, i.e., for the half- 
filled band. 159 However, in k-(BEDT-TTF) 2 X, a differ- 
ence from the cuprates lies in that the band filling (of the 
dimer model) remains at half filling even when metallized 
or superconducting upon increasing the (chemical) pres- 
sure. In this context, Baskaran pointed out the possibility 
of "self-doping" of carriers at half filling, where an equal 
number (No) of doubly occupied sites and empty sites 
hop in the background of singly occupied sites, which is 
shown to be equivalent to the usual t-J model with 2iV 
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Fig. 8. The FLEX calculation results on the original four 
band lattice. 140 The adopted hopping integrals are those for 
X=Cu(NCS)2- (a)Eigenvalues of the Eliashberg equation for 
d xy and d x 2_ y 2 pairings plotted as functions of temperature. 
(b)Contour plot of the spin susceptibility at T/t^i = 0.01. 

holes. 160 On the other hand, Powell and McKenzie stud- 
ied a model that contains the nearest neighbor (which 
will be denoted as J c here) and the next nearest neigh- 
bor ( Jb) superexchange terms in addition to the Hubbard 
model. 161 In this model, the double occupancy of a site 
is not prohibited, so that the hopping term does not van- 
ish at half filling, thereby circumventing the difficulty in 
the usual t-J model, which is always insulating at half 
filling. Assuming the parameter values J c /t c = 1/3 and 
Jb/Jc — (tb/t c ) 2 < 1, a first order transition from d x i_ y i- 
wave superconductivity to a Mott insulating state was 
shown to occur upon increasing U within the Hartree- 
Fock-Gor'kov approximation. 

3.6 Studies for the Hubbard model on the original four 
band lattice 

The d x 2_ y 2 pairing in the dimer model for t^/t c < l 162 
nevertheless contradicts at least with the experimental 
observation of d xy -\\\te gap for X=Cu(NCS)2 in the ther- 
mal conductivity 144 and STM experiments. 143 Motivated 
by this discrepancy between the theories and the ex- 
periments, Kuroki et al. 140 performed a FLEX study 
on the original four band lattice (Fig. 6 (a)), where they 
adopted the hopping integral values obtained from ex- 
tended Hiickel calculation for X=Cu(NCS)2 given in Ta- 
ble I. 120 It can be seen from these hopping integral values 
that the ratio between the intradimer hopping (tbi) and 
the largest interdimer ones (t^) is about 2, 163 which may 
not be considered as so large. 

From the temperature dependence of the eigenvalues 
of the linearized Eliashberg equation for the two types 
of pairing shown in Fig. 8 (a), they concluded that d xy - 
like pairing dominates over d x 2_ y 2. The origin of this 
result can be found in the spin structure. Namely, the 
spin susceptibility Xi shown in Fig. 8(b) peaks around 
Q ~ (±0.47r, ±0.67r), which is a consequence of a par- 
tial nesting between the open and the closed portions of 
the Fermi surface (Fig. 7(b)). As a result, the gap func- 
tion changes sign between the two portions of the Fermi 
surface, but does not change sign within each portion. 

When the dimerization is strong, the four band model 
approaches the single band model, so that the nesting 
shown in Fig. 7(a) dominates and (icy-like pairing gives 
way to d x 2_ y 2 when tbi is large. 140 In fact, further four 



band analysis has shown that whether d xy or d x 2_ y 2 
dominates depends on (a) the strength of the dimer- 
ization, (b) whether tb/t c [= —tb2/(—t p + t q )\ is close to 
unity (which is a measure for how close the system is to 
an isotropic triangular lattice in the strong dimerization 
limit), and (c) the magnitude of the Fermi surface split- 
ting due to t c t' c (i.e., t p ^ t' p , t q ^ t' q ). 164 Namely, d xy 
is more favored when the dimerization is weaker, tb/t c 
is closer to unity (closer to isotropic triangular lattice), 
and the Fermi surface splitting is larger. These factors 
strongly depend on the anions. For X=Cu[N(CN)2]Br, 
the dimerization is stronger {tb\/tb2 ~ 2.7 as compared 
to ^2.0 for Cu(NCS)2), h/t c ~ 0.7 is away from unity 
(compared to ~0.8 for Cu(NCS)2, see Table I), and there 
is no Fermi surface splitting because t c = t' c . In this case, 
d x 2_ y 2-w&ve is expected to dominate over d xy at least 
within the FLEX approach. 

Although the results of the four band approach seem 
to be consistent with the d Ky -wave pairing observed 
in the thermal conductivity 144 and the STM experi- 
ments 143 for k-(BEDT-TTF) 2 Cu(NCS) 2 , there remains 
an issue concerning the T c . Namely, T c exceeding 10K 
for X=Cu(NCS)2 is almost the highest among charge- 
transfer-type organic superconductors. On the other 
hand, the eigenvalue A in Fig. 8 remains small in the tem- 
perature range studied in the FLEX study. Although the 
low temperature regime was not studied in ref.140 due 
to the restriction of the calculation, it is questionable 
whether the eigenvalue actually reaches unity at around 
T ~ 0.004£{,i, which corresponds to the actual T c ~ 10K. 
The eigenvalue becomes larger in the presence of Fermi 
surface splitting (which was not considered in ref.140) 
because the nodes of the d xy gap do not intersect the 
Fermi surface, but the T c problem still exists even in that 
case. 164 In fact, Kondo and Moriya have studied the four 
band model using FLEX, and showed that a realistic T c 
can be obtained only when the dimerization is extremely 
strong (\tf,i/tb2\ > 5). 165 Such a strong dimerization in- 
deed contradicts with the extended Hiickel 119 ' 120 or the 
first principles calculations, 166, 167 and also leads to the 
d x 2_ y 2-w&ve pairing, which is not in agreement at least 
with the thermal conductivity 144 and the STM experi- 
ments 143 for X=Cu(NCS) 2 . 

This problem may be due to one (or more) of the fol- 
lowing possibilities, (i) The FLEX approximation is not 
sufficient for quantitative estimation of T c particularly 
for the K-type BEDT-TTF salts because the nesting of 
the Fermi surface is not good, so that a strong develop- 
ment of a single spin fluctuation mode, which is neces- 
sary to verify RPA-like approaches, is absent especially 
when the dimerization is not so strong. Note that the 
spin fluctuations tend to be weak when the dimeriza- 
tion is weak because the band becomes more closer to 
3/4 filling, where the effect of the on-site repulsion U 
is weaker than for half filling, (ii) In the limit of weak 
dimerization, the band is 3/4-filled and thus far away 
from half filling. In systems away from half filling, terms 
neglected in the FLEX approximation, such as the ver- 
tex corrections, may play an important role, (iii) The 
Hubbard model is oversimplified, and additional terms 
such as the off-site (inter-molecular) repulsions and/or 
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clcctron-phonon interactions are necessary. Note that a 
phonon-mediated attractive interaction that gives large 
contribution around q ~ as considered in section 3.4 
is likely to enhance all types of pairings, so that a con- 
sideration of electron-phonon interaction together with 
the spin fluctuations might result in a superconductiv- 
ity with appropriate T c , maintaining the dominant gap 
symmetry obtained by applying FLEX to the purely re- 
pulsive model. 

Whether the above (i)~(iii) actually provides solution 
to this T c puzzle remains open for future study, 168 but 
some hints have been found in the recent theoretical stud- 
ies on /?'-(BEDT-TTF) 2 ICl 2 , 7 ' 171 ' 172 a superconductor 
which has been found under high pressure by Taniguchi 
et al. 169 /3'-(BEDT-TTF) 2 ICl 2 also has dimcrization of 
molecules, and the T c (^ 14 K) is the highest among the 
charge-transfer type molecular solids and somewhat close 
to the T c of the K-type salts, so that the comparison be- 
tween the two types is intriguing. Using the hopping inte- 
gral values obtained from the first principles calculation 
by Miyazaki and Kino, 170 Kino et al. performed a FLEX 
study on the single band dimer model, and obtained a 
phase diagram in the pressure-temperature space, which 
is similar to the experimental phase diagram including 
the values of T c , although the superconducting phase is 
shifted to a somewhat higher pressure regime. 171 The su- 
perconducting gap in this case has a d xy -like structure, 
reflecting the good nesting of the Fermi surface. Later, 
Nakano and Kuroki 172 performed a FLEX study on the 
two-band lattice with finite dimerization. This was mo- 
tivated by the fact that the ratio between the intradimcr 
hopping integral and the largest interdimer one was esti- 
mated to be about ~ 2, 170 which is a situation similar to 
the case of k-(BEDT-TTF) 2 Cu(NCS) 2 , 120 so that adopt- 
ing the dimer model might be questionable. Nevertheless, 
the gap structure, T c , and the phase diagram obtained 
in the two-band approach is found to be quite similar 
to the one obtained in the dimer model approach. 172 
This suggests that the possibility (ii) mentioned above 
is not, or at least not always, the case. The fact that the 
Fermi surface of /3'-(BEDT-TTF) 2 ICl 2 is well nested may 
have some relevance to the point (i) mentioned above, as 
pointed out in ref.172. 

3. 7 Tests for the Pairing Symmetry Candidates 

Here we survey some theoretical proposals for further 
experimental tests on the pairing symmetry. 1 Powell 
and McKenzie 174 showed that the T c dependence on 
the level of structural disorder 175 ' 176 can be explained 
by assuming d-wave pairing and using the formula (3). 
Since the same applies to s-wave pairing in the pres- 
ence of magnetic impurities, and there is indeed a possi- 
bility that disorder leads to a formation of local mag- 
netic moments due to the proximity to antiferromag- 
netism, the two possibilities (d+non-magnctic defects or 
s+magnetic ones) cannot be distinguished within these 
experiments alone. They further proposed probing the 
presence/absence of magnetic impurities by, e.g., a /LtSR 
experiment. 

Li showed that d xy , d x 2_ y 2, and s-wave superconduct- 
ing states all exhibit different spin structures, so that 



they can be distinguished by looking at the spin sus- 
ceptibility below T c . 177 Namely, for s-wave, a full gap 
opens on the Fermi surface so that the spin susceptibil- 
ity is strongly suppressed. As for the pairings with gap 
nodes, since d xy and d x 2_ y 2-wave pairings open up a gap 
at different portions of the Fermi surface, different types 
of nesting take place (as mentioned in section 3.6, two 
types of nesting are possible for the Fermi surface of the 
k salts), resulting in a different spin structure. 

Tanuma et al. proposed a magnetotunneling spec- 
troscopy for determining the pairing symmetry. 178 
Namely, as mentioned in section 3.3, the density of states 
of a superconductor having nodes in the gap oscillates by 
rotating the direction of the magnetic field. 145 This oscil- 
lation can be detected by looking at the surface density 
of states in the tunneling spectroscopy. The phase of the 
oscillation is different between d x 2_ y 2 and d xy pairings, 
so that the pairing symmetry can be determined as in 
the case of the thermal conductivity measurement. 144 

As in the case of TMTSF salts, it is worth mention- 
ing that the pairing symmetry need not be the same for 
different anions since here again, the competition among 
the pairing symmetries may be close. As mentioned in 
section 3.6, the competition between d xy and d x i_ y i can 
be affected by various factors, such as the splitting of the 
Fermi surface and the strength of the dimcrization, which 
depend on the anions. Speaking of anion dependence, it 
is also important to notice that d xy -wave pairing in k- 
(BEDT-TTF) 2 Cu(NCS) 2 may look like s-wave pairing 
from those experiments that detect only the amplitude 
of the gap because the nodes of the gap do not inter- 
sect the Fermi surface due to the splitting of the Fermi 
surface (see Fig. 7(b)) just like d-wave and /-wave in the 
anion ordered (TMTSF) 2 C10 4 . 

4. Concluding Remarks 

In this paper, a review has been given on the theoret- 
ical studies concerning the pairing symmetry competi- 
tion in (TMTSF) 2 X and k-(BEDT-TTF) 2 X. Some of the 
pairing symmetry candidates and their possible micro- 
scopic origins have been discussed. Existing experimental 
results concerning the pairing symmetry as well as some 
proposals for further tests have also been surveyed. Close 
competition among different pairing symmetries makes 
the problem of theoretically pinning down the symme- 
try difficult, and at the same time, very intriguing. One 
of the origin of this close competition is the peculiarity 
of the Fermi surface, i.e., disconnected Fermi surfaces in 
both (TMTSF) 2 X and k-(BEDT-TTF) 2 X (although in 
different senses), and two types of nesting in k-(BEDT- 
TTF) 2 X. This also implies that several pairing symme- 
tries can be closely competing in the actual materials, 
so that the symmetry might be different for different an- 
ions, or even for the same anion under different envi- 
ronment (magnetic field, pressure, etc.). Although a lot 
has been understood concerning the superconductivity 
in these materials, further theoretical and experimental 
studies are required for a more clear understanding of 
the nature of the superconducting state. 
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